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Abstract 

Previous  investigations  concerning  the  scattering  of  elastic 
waves  from  solid  spherical  inclusions  have  furnished  expressions 
for  the  scattering  cross  sections  which,  upon  numerical  evaluation 
exhibited  resonance-like  features  as  a  function  of  frequency.  In 
the  present  work,  we  study  these  resonances  in  a  fashion  suggested 
by  the  resonance  theory  of  acoustic  scattering  due  to  Flax, 
Dragonette  and  Uberall  £j.  Acoust.  Soc.  Am.  63^,  723  (1978)  J  .  The 
resonances  of  the  solid  inclusions,  exemplified  by  iron  or  lucite 
spheres  imbedded  in  an  aluminum  matrix,  are  found  numerically  in 
the  individual  normal-mode  scattering  amplitudes,  and  are  > 

interpreted  in  terms  of  phase-matched  circumferential  waves. 
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Dispersion  curves  for  the  phase  velocities  of  the  latter  are 
obtained,  exhibiting  two  families  of  waves  of  different  type. 
Finally,  the  connection  of  these  waves  with  the  Stonely  waves 
on  the  boundary  between  two  flat  half-spaces  is  noted  in  the 
high-frequency  limit. 

INTRODUCTION 

The  scattering  of  plane  compressional  (P)  or  shear  (S) 

waves  by  cavities  or  inclusions  in  an  unbounded  elastic  solid 
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medium  has  been  investigated  by  a  number  of  authors  .  For 
example,  the  scattering  from  cylindrical  obstacles  has  been 
studied  for  the  case  of  empty  1  or  fluid-filled^-^  cavities, 
for  both  incident  P  and  S  waves.  Scattering  from  ellipsoidal 
solid  inclusions  has  been  analyzed  for  P  and  S  waves^.  Spheri¬ 
cal  obstacles  have  been  the  object  of  numerous  investigations: 
e.g.,  for  empty  spherical  cavities,  the  scattering  of  P  waves 
has  been  dealt  with  for  the  steady-state^  and  the  transient 
case.  Scattering  from  fluid-filled  spherical  cavities  was 
considered  for  incident  P  waves®”^  as  well  as  S  waves 
An  extensive  review  of  the  literature  on  elastic-wave  scattering 
from  cavities  and  inclusions  is  given  in  Reference  10. 

For  the  case  of  scattering  from  solid  elastic  spherical 
inclusions,  the  literature  contains  numerical  studies  of 
scattering  cross  sections  for  incident  P  waves^  and  S  waves*4, 
f  .  •;  *  *’  V!*>V  ,  |'  l  •*'"  rff***#  •*  <V 
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based  on  the  corresponding  preceding  theoretical  work1'*'11. 

i 

When  plotting  the  cross  sections  against  a  dimensionless 
frequency  variable,  the  authors  of  these  studies  noticed  in 
many  cases  a  "peaking  and  fluctuating"  or  "oscillatory 
behavior"  of  the  plots,  "suggesting  that  some  type  of  resonance 
process  is  in  effect".  No  detailed  analysis  of  the  resonance 
phenomena  had  been  attempted  at  that  time,  however. 

The  present  work  devotes  itself  to  such  an  analysis,  using 
the  recent  theory  of  acoustic  resonance  scattering  due  to  Flax, 
Dragonette  and  Uberall-1-0.  This  analysis  is  carried  out  for  the 
case  of  plane  compressional  waves  incident  on  spherical  elastic 
inclusions  imbedded  in  an  unbounded  elastic  solid  medium. 

(Note  that  a  resonance  study  of  scattering  has  been  performed 
for  cylindrical3,4  and  spherical  fluid-filled  cavities9,10,12'1'7 
A  similar  analysis  is  now  being  carried  out  by  us  for  shear 
waves  incident  on  spherical  elastic  inclusions.  Our  resonance 
theory  separates  the  individual  ("partial  wave")  mode  contri¬ 
butions  into  a  non-resonant  ‘geometrical"  background  term,  and 
into  a  resonant  term  due  to  the  vibrations  of  the  inclusion. 

The  resonant  term  is  small  off  resonance,  but  is  significant 
and  interferes  with  the  background  at  and  near  the  resonance. 

The  position  in  frequency  of  the  resonances  can  be  determined 
by  solving  a  corresponding  eigenvalue  equation,  which  for  the 
case  of  non-absorbing  media  is  a  real  equation  with  real 


solutions . 


We  analyze  numerically  the  cases  where  the  scatterer  is  (a) 
a  lucite  sphere,  and  (b)  an  iron  sphere,  both  imbedded  in  an 
aluminum  matrix.  These  materials  were  chosen  in  order  to 


provide  a  contrast  in  the  corresponding  background  terms 
^approximately  given  by  that  of  (a)  an  empty  cavity,  or  (b)  a 
rigid  inclusion}.  The  resonances  are  obtained  numerically, 
and  are  interpreted  in  terms  of  phase-matched  circumferential 
waves.  Dispersion  curves  for  the  latter  are  obtained, 
exhibiting  two  families  of  waves  of  different  type.  The 
connection  of  these  surface  waves  with  the  Stoneley  waves  on 
the  boundary  between  two  flat  half-espaces  is  noted  in  the 
high-frequency  limit. 


I.  THEORY 

The  center  of  an  elastic  sphere  of  radius  a  is  taken  as 
the  origin  of  a  spherical  coordinate  system  (r,-$,cp)‘.  The 
longitudinal  and  shear  speeds,  and  the  density,  are  designated 
by  cL,  cs  and  ^  ,  respectively,  for  the  elastic  matrix  in 

which  the  obstacle  is  imbedded,  and  by  v^,  Vg  and  g ^  for  the 
spherical  obstacle  itself. 

The  displacement  vector  u  in  either  elastic  material  can 
be  expressed  in  terms  of  two  disjJacement  potentials  J5"  and  , 
i.  e 


u  =  -  V  <£>  +Vx^ 


(1) 


where  S2  is  the  scalar  potential  and  the  vector  potential. 

Each  potential  satisfies  a  scalar  wave  equation.  For  the  incident 
wave,  we  take  as  its  solution  a  plane  harmonic  longitudinal  wave 
propagating  in  the  direction  of  the  z  -  axis,  with  wave  number 
kL  and  circular  frequency  CJ.  This  wave  is  given  by 
2?  *  -  expi(to.t  -  k  lZ-)  which  can  be  expanded  in  spherical 

A. 

coordinates  as 
00 

=  ^"(2n  +1)  (-‘i)njn^Lr^Pn^cos'^  )  *  '  (2) 

sn-o 

As  a  result  of  the  scattering  and  mode  conversion  processes, 

there  are  both  longitudinal  and  shear-type  scattered  waves,  i.e. 
ext 

<2n  +  1)  (  ~i)n  An&n(2*  (kLr)  Pn(cos^>)  (3a) 

'Yl'O 

(longitudinal) ,  and 
o o 

(2n  +  1)  <  -i)nBnAn(2>  (^rj^tcos^)  (3b) 

/»v=o 

(shear),  where  by  jf  s  we  designate  the  <j p—  component  of 
y  g  which  is  the  only  nonvanishing  one.  The  wave  numbers  are 
here  kL  =&)/ cL  and  k.s  =  C0/cs;  hn^)  s  Hankel  I 

functions  of  the  second  kind,  and  Pn,  Pn^  are  the  ordinary 
and  associated  Legendre  polynomials.  The  undetermined  coefficients 
An  and  Bn  have  to  be  determined  from  the  boundary  conditions  at 
r  =  a. 
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Inside  the  scatterer,  the  refracted  waves  (labeled  by  a 
subscript  q)  are  represented  by 

DO 

S  =  £>n  +  l)<-i>nCnjn<tIr)Pn  (  COS  $)  (4a) 

and 

oo 

')}?  =  (2n  +  (cos^t  )  ,  (4b) 

*  '*X=-0 


where  the  wave  numbers  are  qL  =  C<yvL  and  qs  =  &VVg. 

There  are  four  separate  boundary  conditions  that  must  be 
satisfied  by  the  solution  to  this  problem.  These  conditions 
are  prescribed  at  the  interface  r  =  a  between  the  sphere  and  the 
matrix;  they  are:  continuity  of  (i)  normal  (radial)  and  (ii) 
tangential  displacements,  and  continuity  of  (iii)  normal  and 
(iv)  tangential  stresses.  These  are  sufficient  to  specify  the 
four  sets  of  unknown  coefficients  An,  Bn,  Cn  and  Dn.  Details 
of  setting  up  these  boundary  conditions  are  given  in  the 
literature^S  and  we  shall  here  simply  state  the  results. 

The  coefficients  of  the  scattered  waves  are  given  as  the 
quotients  of  two  fourth-order  determinants  of  the  form 


all  ai2  a14  a15 
a21  a22  a24  a25 

a31  a32  a34  a35 


(5a) 


J. 

A 


a41  a42  a44  a45 
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II 


and 


B 


n 


A  - 


all 

a12 

a13 

a15 

1 

a21 

a22 

a23 

a25 

A 

a31 

a32 

a33 

a35 

a41 

a42 

a43 

a45 

denominator  is 

all 

a12 

a13 

a14 

a21 

a22 

a23 

a24 

a31 

a32 

a33 

a34 

a41 

a42 

a43 

a44 

(5b) 


(5c) 


The  twenty  different  elements  contained  therein  are  found  to  be 


all  = 


l 


?  i-  *-1 


[  A  +  ^  fa"  (6>i- )] 


l12 


*13 


=  _2f 


(6a) 


*14 


oo-fx+l) 


(from  the  continuity  of  normal  stresses  t  ) , 

i-r- 

c*^=  °L.tf2  (®u) 

<321.  =  (Qg) 

^  =  =zl£:»'  czu)  < 

a2i  =  00(71+1) -£JX) 

Ol5r  =  Jx  (iL) 

(from  the  continuity  of  normal  displacements -to  ), 

*7 

^  ^  fQL  W-- (QJ]  - 

<'*32-  =  ^  ^  /^j)  +  (nta)(n-i)^HfQj  )^ 

^,J5  =  ^j^'rzO-^rzj] 

=  2*££)“(ZS) i-  (*-,+z )fH-<t)^V-2^) 


(6c) 


-a* --2[2^teL)-^CZL)] 

(from  the  continuity  of  tangential  stresses  ,  and 

fQ$)  ^\a  fQ$) 

~  r^L ) 

/0.4£r  *  ^ 


(from  the  continuity  of  tangential  displacements  ) . 


Here 


,  kLa,Zg  kL a,  QL-^L^y  Qy- 


C  *“  /■*-£.  / ^%  )  y  ~  whereA/yU- 

the  Lame  constants  of  the  matrix  and  A  u.  those  of  the 

4  ’r<i. 


sphere,  so  that 
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-cL=  0+ir)/?3''*- 

sCf,  =  <7al 

vl  =  [  ( A  J  +  2  H  ^  (7b> 

'u?=  (H/?+V*u  • 

This  determines  the  coefficients  of  the  scattered  wave  in  terms 
of  the  incident  amplitude. 

The  quantity  of  physical  interest  is  here  the  differential 
scattering  cross  section  as  obtained,  e.g.,  by  Ying  and  Truell^^. 
The  time  rate  at  which  energy  is  carried  away  by  the 

scattered  wave  across  a  sphere  of  radius  R  OO  is 

7-"  =  J  (Trr  a*  +  "2"^ <8a> 

with  real  displacements  u^  and  stresses  T^.  .  With  complex 
fields,  and  going  over  to  spherical  components,  this  gives 

*  =  jr£e  (yrr  u*  +  )TZ1dSL  <sb) 

since  there  is  no  Cp  -  dependence  present;  Using  the  well-known 
expressions  for  displacements  and  stresses  in  terms  of  the 
potentials  jjT  and  ,  e.g.  Eq.  (1),  and  inserting  the  far-field 


expressions  of  Eqs.  (3a)  and  (3b),  we  find 


T 


T 


-  11  - 


dr  -  i  usL  fu>*  { i  f  fjV)  r+ 1  £  *>>n i9a> 

L  rt=o  -vv"  i 


where 

a" 4/a-  <^)  —  Cx-w+'il^L  A^  r<os-^j 

a~1/a  £?*($)  =  (2^  +  1)^  "B^  T^fros^). 


' 

(9b) 


(9c) 


The  total  cross  section  6*  is  obtained  as  the  integrated 
flux  F,  normalized  to  the  flux  (i/2)^o>2kL  of  the  incident 
wave,  so  that 


C r 

7caa 


0v~  o  L- 


2. 

+  'Yl('Yl  +  'l) 


(IQ) 


Note  that  for  backscattering  (S=k) 
to  the  scattering  amplitude  from  the 

n 

(ftr)  since  Pn1  (-1)  «  0. 


there  is  no  contribution 
Outgoing  shear  waves 
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II.  RESONANCE  THEORY 


The  resonance  theory  of  acoustic  and  elastic  wave  scattering 
has  been  described  in  detail  in  references  16  and  9, 
respectively,  and  shall  only  be  sketched  in  the  present  context. 
Using  the  asymptotic  forms  of  the  Hankel  functions,  one  obtains 
the  far-field  expressions  of  Eqs.  (3)  as 

&  ~  (i/r)  ( rn  e ( ^  fT  ■  <ll*> 

^  ~  H/y)  (  aw 

OA. 

together  with  Eqs.  (9b, c) .  TWo  S-matrix  elements  are 
introduced  by 


and  we 


H  ?T> 

0-V 

(12a) 

S™  - 

2B 

• 

ov 

(12b) 

find 

".I 

A 

i 

Q  Id  ^21 

(13) 

' 

A 

A 

A 

■^•41  ^42. 

i 


-  13  - 


where  equals  -  but  with  hn(2)  replaced  by  hn(1) 

It  is  advantageous  to  consider  the  limit  of  scattering 
from  a  "soft"  (i.e.  evacuated)  spherical  cavity.  This  is 
obtained  by  ^ -*  0  while  vL  and  vg  remain  constant,  so  that 
remains  finite  and  0.  (The  "rigid"  limit  would  be 
obtained  by  )  . 

Ohe  corresponding  "soft"  S  matrix  element  is  thus 


(s)TP 


c*.n 

<»vi 

In  the  resonance  theory9'16,  this  expression  is  factored  out  of 
PP 

Sn  and  we  can  write 


< ^  ^  (s)TT  6?^  _  (s)TP  +•  *  1y*  G 

i?tC,  t  *•  q*. 


(15a) 


with 


d  = 


^Xi. 

A 

A 

oi23 

a;i 

a3l 

A 

oHZ 

(15b) 
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and  Gn  a  similar  expression  without  carets.  The  resonances  in 
the  scattering  amplitude  originate  from  the  condition  (repre¬ 
senting  a  real  eigenvalue  equation) : 


Re  Gn(2*L)=  0, 


(16) 


which  determines  the  eigenfrequencies  ZL=  ZLnr  labeled  by 
r  «  1,2,3...  .  A  linear  expansion  of  G  about  then  brings 

Snpp  into  an  explicit  resonance  form: 


5"  »  sir”  E 


OO 


H- 


>ir 


=i 


with  a  resonance  width 


and  with 


^  -  4 /-Re Q^'CZ^) 


(17a) 


=  2.1-mG^  (2r)/Re  6^  ( <-> 


(17c) 


If  Eq.  (14)  is  written  as 


TP 


(*> 


2* 

with  complex, 

An  =  |  (Sn*P-  1) 


then  the  quantity 


(18a) 


(18b) 


-  15  - 


which  determines  the  scattered  amplitude  of  Eq.  (3a) ,  becomes 


r«i  rw  /2. 


4-  2-  >c  G 


SiVl 


]  . 


(19a) 


where 


=  (  fix  -  ). 


(19bl 


This  shows  that  the  scattering  amplitude  consists  of  a 
non-resonant  background  term  £the  second  term  in  brackets  in 
Eq.  (19a)^  which  is  due  to  the  scattering  of  P  waves  from  a 
soft  sphere,  with  a  series  of  superimposed  resonances  (first 
term  in  brackets)  which  will  interfere  with  the  background 
amplitude.  A  corresponding  expression  may  be  derived  for  Bn 
of  Eqs.  (3b)  and  (12a),  i.e.  the  mode-converted  scattering 
amplitude: 


-  16  - 


£  £ 
2.  ^ 


i  e2.*r. 


if 


£  VC^r/i 


ij 


(19c) 


Note  that  the  resonance  frequencies  are  here  the  same  as  in 
Bq.  (19a)  since  the  resonances  are  due  to  the  vibrations  of  the 
same  inclusion  in  both  cases,  independent  of  the  resonance  • 
mechanism. 

Ihe  expressions  of  Eqs.  (19a),  (19c)  hold  in  the  framework 
of  the  resonance  approximation,  valid  if  the  width  of  the 
resonances  is  samll  compared  to  their  spacing.  However,  the 
actual  numerical  calculations  to  be  discussed  below  have  all 
been  carried  out  using  the  exact  expressions  Eqs.  (5) . 


III.  NUMERICAL  RESULTS 


Scattering  amplitudes  have  be  evaluated  numerically  for  the 
cases  of  (a)  a  lucite  sphere,  and  (b)  an  iron  sphere,  both 
imbedded  in  an  aluminum  matrix.  The  material  parameters  of 
these  substances  (assumed  non-absorptive)  are  listed  in  Table  I. 

In  Pig.  la,  we  plot  for  case  (a),  lucite  in  aluminum,  the 
amplitude  moduli  .«x | /I^  (co^iS)of  Eq.  (9b)  kLd  = 


of  aluminum,  for  the  partial  .  waves  n  =  0,  1  and  2.  Similarly, 
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Pig.  lb  shows  the  quantity  a~1//2  I  ^nPS  I  /Pn^  (cosiMof 
Eq.  (9c).  i.e.  the  mode-converted  scattering  amplitude.  In 
both  cases,  the  interference  pattern  of  the  resonances  of 


the  lucite  sphere  with  a  non-resonant  background  is  clearly 
visible.  Due  to  the  much  smaller  impedances  of  the  lucite 
sphere  (  *  P-fc  ^  )  as  compared  to  those  of  the  aluminum 

matrix  (  ^  ) ,  it  is  expected  that  the  background 

resembles  closely  that  of  scattering  from  an  empty  cavity  in 
aluminum.  We  therefore  plot  in  Figs.  2a  and  2b  the  quantities 

a  'lA  If”'  C”  l/T*.  and  «.-*■  I/J*  -  I  /*?  , 

respectively,  where  the  "soft"  background  amplitudes 


'  obtained  from  Eqs.  (9b),  (9c)  in  the 
limit  0,  were  subtracted  out  from  the  total  partial-wave 

amplitudes.  Indeed,  the  pure  resonances  of  the  cavity  contents 
appear  to  have  been  made  visible  in  this  way.  The  resonance 
frequencies  are  the  same  in  both  cases,  but  it  is  seen  that 


rnT-S 


,  obtained  from  Eqs.  (9b),  (9c)  in  the 


two  families  of  resonances  are  present,  and  only  one  of  these 
is  being  excited  in  P  ->  P  scattering  while  both  families  are 
excited  in  P->  S  scattering. 

Similar  numerical  results  are  shown  for  case  (b) ,  iron  in 
aluminum,  where  the  n  =  0,  1  and  2  partial  wave  scattering 
amplitude  moduli  are  shown  for  P  P  scattering  in  Fig.  3a,  and 

for  P  S  scattering  in  Fig.  3b.  Due  to  the  density  of  the 
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inclusion  being  larger  than-  that  of  the  matrix,  the  background 
is  now  expected  to  correspond  to  scattering  from  a  rigid  sphere 
(^>/p^  O)  .  In  fact  if  this  background  is  bubtracted  from 

^  or*  and  the  modulus  is  taken  subsequently.  Figs.  4a 

(P  P  scattering)  and  4b  (P  -*  S  scattering)  indicate  that  the 
pare  resonance  amplitudes  are  indeed  obtained  in  this  way. 

(The  resonance  formulation  with  a  rigid  background  is  similar  to 
that  for  the  soft  background  given  in  Section  II,  see,  e.g. 
references  9  and  16) . 

From  Figs.  2  and  4,  the  resonance  frequencies  for  cases  (a) 

and  (b)  may  be  read  off  and  tabulated.  They  may  serve  to 

obtain  the  dispersion  curves  of  circumferential  waves  propagating 

over  the  surface  of  the  inclusion,  which  have  been  shown  to  be 

generated  during  the  scattering  process,  and  to  provide  an 

interpretation  for  the  existence  of  the  resonances10'19. 

The  resonance  frequencies  2Lnr  in  the  resonance  expression 

for  the  S  -  matrix,  Eq.  (17a),  are  functions  of  the  mode  number 

n  which  we  now  consider  a  continuous  variable.  We  may  expand 
— _  ■nr' 

in  the  n-  variable  about  the  value  n  =  nr,  chosen  so  that 
(the  incident  frequency)  : 


Z”  =  ZLn,’r+  K-,v)  (d.'E^/d.rv) 


(20) 
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X  exp  l[2  -  7r  >V  /2  +  1  )7frir  +  £ I  )  .  (22 

This  describes  a  manifold  (labeled  by  r)  of  attenuated 
circumferential  "tidal"  waves  which  engulf  the  spherical 
inclusion  along  all  meridians  (m  =  o) ,  with  additional  components 
(m  y  O  ) that  have  previously  circumnavigated  the  inclusion  m' 
times.  The  phase  velocity  of  the  rth  surface  wave  is  given  by 


When  the  Regge  pole  of  Eq.  (21d) ,  which  with  varying  frequency 
moves  through  the  complex  n  -  plane  along  a  "Regge  trajectory", 
passes  by  the  integer  n  so  that  nr ->  n,  the  wavelength  of  a 
surface  wave, 

A^n=L)  =  2irot/r-nr  ). 

becomes  2xoi/(rc  +  i)  so  that  n  +  ~  wavelengths  span  the 
circumference  of  the  spherical  inclusion.  Since  as  seen  from 
Eq. (22) ,  the  circumferential  wave  suffers  a  phase  loss  of  Tt/2 
(i.e.,  a  quarter  wavelength)  every  time  it  passes  through  one  of 
the  convergence  points  *$  =  0,7T*  at  the  north  and  south  pole  of  the 
spherical  inclusion,  this  condition  leads  to  a  perfect  phase 
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match  of  the  wave  after  each  circumnavigation.  The  poles  in  the 
scattering  amplitude,  Eqs.  (17a)  or  (21a),  therefore  describe  a 
"resonant  buildup"  of  the  amplitude  of  the  circumferential  wave 
in  the  course  of  its  repeated  circumnavigations  due  to  phase 
matching,  which  becomes  physically  manifest  at  the  integer  mode 
amplitudes,  leading  to  the  resonances  shwon,  e.g.,  in  Figs. 2  or  4. 
The  families  of  resonances  which  in  these  figures  are  manifest 
in  all  partial-wave  amplitudes  and  which  recur  at  successively 
higher  frequencies  in  successive  partial  waves,  each  thus 
correspond  to,  and  represent,  one  single  circumferential  wave. 

In  a  sense,  this  physical  pictOre  (which  applies  generally  to 
any  type  of  smooth  spherical  scatterer) ,  explains  the 
existence  of  the  scatterer' s  resonances  in  terms  of  a  phase 
matching  of  the  circumferential  waves  which  its  surface  may 
support. 

Using  the  resonance  frequencies  which  appear  in  Figs.  2  and  4, 
(as  well  as  in  corresponding  graphs1  for  higher  partial 
waves) ,  we  may  obtain  dispersion  curves  of  the  surface  waves 
from  Eq.  (23)  .  These  are  shown  in  Fig.  5  (for  lucite  in  aluminum) 
and  Fig.  6  (for  iron  in  aluminum),  where  the  quantity  cr/cL  is 
plotted  vs.  ZL  for  all  the  surface  waves  labeled  by  r.  From 
these  dispersion  curves,  we  see  that  the  surface  waves  on  an 
elastic  spherical  inclusion  imbedded  in  an  elastic  medium  fall 
into  two  classes  (represented  in  Figs.  5  and  6  by  solid  and 
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dashed  curves,  respectively),  which  presumably  correspond  to 
longitudinal  and  shear  type  surface  waves.  The  distinction 
between  these  two  classes  of  curves,  when  producing  the  plots, 
has  been  greatly  facilitated  by  the  fact  that  in  P->  P 
scattering  (Figs.  2a  and  4a)  only  one  cla  i  of  surface  waves 
seems  to  have  been  excited,  while  in  P*>  S  scattering 
(Figs.  2b  and  4b),  both  are  excited.  For  small  values  of  Z L, 
however,  the  distinction  in  the  excitation  mechanism,  and  hence 
in  the  dispersion  curves,  is  no  longer  so  clear-cut. 

At  the  end,  we  would  like  to  comment  on  the  high-frequency 
(or  large-a)  limit  of  the  dispersion  curves,  .  In  this 

limit,  the  dispersion  curves  approach  horizontal  asymptotic 
values,  while  physically,  the  boundary  of  the  inclusion  tends 
towards  the  flat  boundary  between  two  infinite  elastic  half-spaces. 

In  that  latter  case,  the  speed  of  surface  waves  on  the  boundary, 
c,  satisfies  the  Stoneley  equation^0 

-+  2."R  'c1  -  f  +  Pf  )  ^ 

■*  "R5-  (MN  -1)  (Mt  -1  )  =  0 
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where 

N  =  (i.  -  'C2-/^L.MVi 

M„  =  (1  - 't2  I^Y11- 

*1 

N  =  ( 1  - 'CI-/'cs'-')1/‘-  < 

=  (  1  -  ^  /V>^')W3' 

"R  =  2  (? 

We  have  solved  this  equation  numerically,  and  have  found  a 

Stoneley  root  c/c^  =  0.504  for  the  case  of  iron  in  aluminum; 

no  root  was  found  for  the  case  of  lucite  in  aluminum.  (It  is 
20 

well-known  that  Stoneley  wave  solutions  do  not  always  exist 
on  the  flat  boundary  between  two  elastic  half-spaces) . 

In  Fig.  5,  we  have  entered  two  arrows  on  the  right-hand  side 
representing  the  values  vl/cl  =  0.4080  and  Vg/c^  =  0.2010. 

While  many  of  the  dispersion  curves  on  this  graph  fall  below 
vl/cl  for  large  enough  ZL  the  value  vs/cL  is  found  to  represent 
a  comfortable  lower  limit  for  even  the  lowest  dispersion  curve, 
up  to  the  largest  ZL~  value  (ZL  =  25)  at  which  that  dispersion 
curve  was  obtained;  but  this  value  was  not  yet  large  enough  to 
decide  whether  vs/cL  actually  constitutes  an  asymptote  for  the 
dispersion  curves. 
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In  Fig.  6,  we  have  entered  three  arrows  representing  the 
values  (from  top)  vl/cl  =  0.8815,  c/cL  =  0.504  (the  Stoneley 
wave  solution),  and  vs/cl=  0.4720.  While  either  of  the  two 
latter  values  seems  to  represent  an  asymptote  to  the  dispersion 
curves  on  this  figure,  a  plot  of  the  lowest  dispersion  curve 
up  to  ZL  =  35  reveals  that  at  ZL  =  24,  this  curve  falls  below 
the  value  of  c/c^  ,  but  that  the  value  of  vs/c l  represents  a 
likely  asymptote  for  the  lowest  dispersion  curve.  Since  in 

this  case  of  iron  in  aluminum,  the  higher  dispersion  curves  lie 

< 

considerably  above  the  lowest  one,  no  statements  could  be  made 
about  their  asymptotics,  except  that  several  of  them  definitely 
fall  below  the  value  of  vL/cLas  Z^  increases. 

IV.  CONCLUS ION 

This  work  represents  the  first  systematic  study  of  resonance 

effects  in  the  scattering  of  elastic  (P)  waves  from  solid 

spherical  inclusions.  The  resonance  approach  is  patterned  after 

16 

our  previous  resonance  studies  in  acoustic  scattering  ,  and  in 
elastic-wave  scattering  from  fluid-filled  cavi ties^ • ^ .  Both  in 
that  latter  case,  and  in  the  present  case,  the  resonarc  es  are 
caused  by  the  vibration  of  the  substance  filling  the  cavity,  and 
they  can  be  shown  to  originate  from  the  phase-matching,  or 
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resonant  reinforcement,  of  repeatedly  circumnavigating  surface 
waves.  But  while  in  the  case  of  a  fluid-filled  cavity  only 

Q  IQ 

one  type  of  surface  wave  has  been  found  '  ,  the  present  case 

of  a  solid  inclusion  is  endowed  with  two  different  classes  of 
surface  waves,  as  evident  from  their  dispersion  curves  which 
were  obtained  here.  It  is  interesting  to  note  that  P->  P 
scattering  (without  mode  conversion)  only  sets  up  one  of  these 
two  classes  of  surface  waves,  while  P  -»  S  scattering  (with  ' 
mode  conversion)  excites  both  classes.  The  asymptotic  behavior 
for  high  frequencies  of  the  dispersion  curves  has  been  compared 
with  the  Stoneley  wave  limits  for  two  half-spaces  in  contact, 
and  it  was  found,  at  least  for  the  case  of  an  iron  inclusion  in 
aluminum,  that  the  Stoneley  wave  speed  does  not  represent  an 
asymptote  of  the  lowest  dispersion  curve  while  the  shear  speed 
of  iron  probably  does.  In  the  lucite-aluminum  case,  again  the 
shear  speed  of  Lucite  probably  represents  the  asymptote  of  some 
of  the  dispersion  curves. 


TABLE  I 


Material 


Aluminum 


Material  Parameters  Used 


density 

xlO~^kq/m^ 


Longitudinal  speed 
x  10  m/sec 


2.7 


6.568 


Shear  speed 
—3 

x  10  m/sec 

3.149 
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FIGURE  CAPTIONS 


Fig.  1.  Plot  of  Ma-,/‘-\f”\/Inr«.,V.  and  (b) 

I  (cos-»0  vs,  2l  for  the 

partial  waves  n  =  o,  1,  and  2,  case  of  a  Lucite 
sphere  in  aluminum. 

Fig.  2.  Plot  of  (a)  ~ I  /  o5  $•) 

and  (O  XX'*-17-  _  ^')T5  1  /P^(cost)  Vf.  "21 ^ 

(where  ^  is  the  empty-cavity  scattering  amplitude) 
for  the  partial  waves  n  =  o,  1,  and  2,  case  of  a 
lucite  sphere  in  aluminum. 

Fig.  3.  Plot  of  (a)  If+J*  I  /P*.  (cot£)/ 

and  (b)  AX~t/z  l |  /  P-*^  (cos$)  VS.  PrL 

for  the  partial  waves  n  =  o,  1,  and  2,  case  of  an  iron 

sphere  in  aluminum. 

Fig.  4.  Plot  of  (a)  ~  )T?  I  /T^(Cos^) 

and  (b)  ^-‘d'/1-  ***  I  /P^1  fcoS^) 

vs.  Z^  (where  is  the  rigid-inclusion  scattering 

amplitude)  for  the  partial  waves  n  =  O,  1,  and  2,  case 
of  an  iron  sphere  in  aluminum 

Fig.  5.  Dispersion  curves  for  the  two  families  of  surface  waves 
on  a  lucite  sphere  in  aluminum.  Arrows  on  the  right 

'  represent  the  values  vL/cL=  0.4080  and  vg/cL  =  0.2010. 

Fig.  6.  Dispersion  curves  for  the  two  families  of  surface 

waves  on  an  iron  sphere  in  alum  in/am.  Arrows  on  the 
right  represent  the  values  Vj/cj,  =  0.8815,  c/cL=  0.504 
(the  Stoneley  wave  spped) ,  and  vs/cL  =  0.4720 


